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We investigate the effects of nonstandard four-fermion neutrino-neutrino interactions on the flavor 
evolution of dense neutrino gases. We find that in the regions where the neutrino-neutrino refractive 
index leads to collective flavor oscillations, the presence of new neutrino interactions can produce 
flavor equilibration in both normal and inverted neutrino mass hierarchy. In realistic supernova 
environments, these effects are significant if the nonstandard neutrino-neutrino interaction strength 
is comparable to the one expected in the standard case, dominating the ordinary matter potential. 
However, very small nonstandard neutrino-neutrino couplings are enough to trigger the usual col- 
lective neutrino flavor transformations in the inverted neutrino mass hierarchy, even if the mixing 
angle vanishes exactly. 
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I. INTRODUCTION 

New neutrino interactions are predicted by several 
extensions of the standard electroweak theory. These 
new interactions can be treated in an effective (low- 
energy) framework by using four-fermion operators, e.g., 
Oq/3 ~ pal^PLV(3\[f^^,PLf] with strength Gap, induc- 
ing either flavor-changing [a ^ (3) or diagonal but flavor 
non-universal neutrino transitions. Nonstandard interac- 
tions (NSI) of neutrinos with charged fermions in mat- 
ter, and their interplay with the oscillation phenomenon, 
have been investigated in many different contexts. An 
incomplete list includes analyses related to the solar neu- 
trino problem ^ 0, 0, 11], to the atmospheric neutrino 
anomal y B, [j , M Bl , to supernova neutrinos [lOl. [ill, [l^. 
El, 0, llSl. Il6|. to primordial neutrinos [l3|, to the pro- 
duction or detection of laboratory neutrinos [3 [l^ , and 
to future long-baseline projects [13, HH [12] ■ Neutrino- 
neutrino interactions are even more difficult to constrain 
experimentally [1^ [13, [1^. In particular, four-fermion 
(left-handed) neutrino-neutrino interactions as large as 
the Standard Model neutral current couplings are viable 
without violating astrophysical or laboratory bounds [26l] 
(see also Ref. 27] for a recent discussion on how to gen- 
erate this type of interactions in a gauge invariant way). 

Dense astrophysical environments like the Early Uni- 
verse, core-collapse supernovae and accretion disks of co- 
alescing neutron stars are among the few systems where 
neutrino-neutrino interactions play a role, via the refrac- 
tive index that can affect the flavor evolution of the sys- 
tem. In particular, it has been shown that flavor oscilla- 
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Since physics beyond the Standard Model can plausibly 
induce at least small deviations from the standard pre- 
dictions, it seems interesting to investigate the impact 
of new neutrino-neutrino interactions on collective flavor 
oscillations in dense neutrino gases. 

We devote our work to this purpose. In Section [ill 
we generalize the equations of motion for a neutrino en- 
semble in the presence of new neutrino-neutrino interac- 
tions. We present our results for the flavor evolution of 
a homogeneous and isotropic gas of neutrinos with de- 
creasing density in Section lllli where we find that NSI 
among neutrinos could completely equilibrate the flavor 
content of the system in both normal and inverted neu- 
trino mass hierarchy. In the case of neutrinos streaming 
off a supernova (SN) core, these new effects are inhib- 
ited by the strong ordinary matter potential unless the 
strength of the NSI is of the same order of magnitude 
as the one predicted by the Standard Model. However, 
also in the case of small values of nonstandard neutrino- 
neutrino couplings, these are enough to trigger collective 
flavor transformations in inverted mass hierarchy for van- 
ishing mixing, where in the standard case no evolution is 
expected. Finally, in Section HVl wc comment our results 
and conclude. 



II. EQUATIONS OF MOTION WITH 
NONSTANDARD NEUTRINO-NEUTRINO 
INTERACTIONS 
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Flavor oscillations of a homogeneous ensemble of neu- 
trinos and antineutrinos are described by an equation of 
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motion (EOM) for each mode p [63| 

= [{nl + v + n^ 



p/>Pp] > (1) 
Pp = [{nl-v-nl),-p^] , (2) 



where [•,•] denotes the commutator. For ultrarelativis- 
tic neutrinos of momentum p, the matrix of vacuum 
oscillation frequencies, expressed in the mass basis, is 
f2p = diag(TOf , m|, m|)/2p, mi=i^2,3 being the neutrino 
masses. At leading order, the only non- vanishing clement 
of the matter potential matrix in the weak interaction ba- 
sis is due to the forward scattering on background elec- 
trons, leading toV = v^Gi^ne diag(l, 0, 0), where Gf is 
the Fermi constant and Ue is the net electron number 
density. Neutrino-neutrino self-interactions introduce an 
additional contribution fip to the refractive energy shift. 
One finds [Hi 
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where dq = d'^q/ {2tt)^ and Vp is the velocity. The factor 
(1 — Vq • Vp) implies "multi- angle effects" for neutrinos 
moving on different trajectories In the Standard 

Model, the dimensionless coupling matrix G is the iden- 
tity matrix. If we consider the possibility of new physics 
beyond the Standard Model, the coupling matrix can as- 
sume a non-universal and/or non-diagonal structure. 

We now focus on a two- flavor system {vctI^x}, where we 
expand all the 2x2 matrices in the EOM in terms of the 
2x2 unit matrix I and the Pauli matrices cr. Explicitly, 
we define [391 
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Here, the overall neutrino (antincutrino) density is given 
by J dp /p — fii, (J dp /p — Uij). For simplicity and with- 



out loss of generality, we here assume that initially only 
Vf. and are present with an excess neutrino density of 
riy^ = (1 + C) "^Pe ; numerical examples, we shall assume 
that the asymmetry parameter between neutrino species 
is ? = 0.25 42]. The vectors Pp and Pp are the neutrino 
and antincutrino polarization vectors. We define the to- 
tal polarization vectors P — JdpPp, P = / dpPp, 
which are initially normalized such that P(0) = (l-|-?)e2 
and P(0) = e^, 62 being the unit vector in the positive 
z-direction. Here, we have chosen our coordinate sys- 
tem in such a way that a polarization vector pointing 
in the positive z-direction represents electron neutrinos, 
whereas an orientation in the negative z-direction cor- 
responds to a combination of muon and tau neutrinos, 
which we denote Vx- We also have o^o — {m\ + m2)/2E, 
the vacuum oscillation frequency is ujp = {ra\ — m\)/2E 
with the energy _E = |p|; L is a unit vector pointing in 
the direction singled-out by the neutrino potential in the 
charged fermion background, and A its normalization; for 
the cases considered here, L = and A = -\/2 Gp^e- The 
unit vector B points in the mass eigenstate direction in 
flavor space, such that B • L = — cos 20, where 9 is the 
vacuum mixing angle. The neutrino-neutrino interaction 
couphngs are given by {50 j g}, with the Standard Model 
case corresponding to Qq = 2, |g| = 0; it is thus the vec- 
tor g which is characteristic of NSI, its third component 
setting possible non-universal couplings (i.e., v^ — Ve dif- 
ferent from Vx — i^x) while its first two components char- 
acterize flavor-violating operators. Note that we have 
assumed that the effective four-neutrino vertex is a good 
description of the new dynamics. This is done having in 
mind an effective field theory correction to the Standard 
Model dynamics induced by an energy scale M 3> M^, 
with |g| - {Mz/Mf. 

With above definitions, the neutrino EOMs assume the 
form 



Pp = (c^pB + AL + nl) X Pp , 



(5) 



where the neutrino-neutrino interaction "Hamiltonian" 

nl is 



^^P = M / dq(l- Vq- Vp) 



50? + g • (Pq - Pq) 
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and we have defined the parameter /i — \/2Gf?t-p which 
normalizes the neutrino-neutrino interaction strength. 
For antineutrinos, the EOMs are the same as for neu- 
trinos with the substitution u-p — > — Wp. It is also trivial 
to check that in the limit of Standard Model-only cou- 



plings (.go = 2, |g| = 0), no direction is singled out in 
the flavor basis and one recovers the standard neutrino- 
neutrino Hamiltonian. 

In our analysis we will assume the single-angle approx- 
imation, in which we substitute the angular structure of 
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neutrino-neutrino interactions in Eq. ([6]) with an effective 
neutrino- neutrino interaction strength /i^ such that: 



Mr 



(7) 



This approximation has been numerically shown to be 
a good description of the flavor evolution in the super- 
nova environment, where a significant asymmetry be- 
tween neutrinos and antineutrinos is expected. We have 
explicitly checked that the single-angle approximation is 
a good description of the flavor evolution in presence of 
nonstandard interactions also for the examples shown in 
this work. In particular, the asymmetry between neutri- 
nos and antineutrinos prevents possible effects of multi- 
angle decoherence in the neutrino ensemble [l^l . Inspired 
by this case, in our numerical examples we use an effec- 
tive neutrino interaction strength [39| 
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where rio = r/(10 km). We consider a spherically sym- 
metric system in which neutrinos are emitted by a sphere 
at i? = 10 km. In this situation, the only spatial variable 
that characterizes the flavor evolution is the radial coor- 
dinate r. Concerning the other input parameters, in our 
numerical examples we will consider a single mode system 
with Lo = 0.3 km~^, corresponding to typical supernova 
neutrino energies for the atmospheric mass squared dif- 
ference. We choose 9 — which, in the standard case, 
would mean that no flavor evolution can occur. We al- 
ways assume go — 2, which essentially means that the 
overall strength of the neutrino-neutrino interactions are 
given by the Standard Model. 

There is an important property which is worth noting: 
If g II L, then the quantity B • D = B • (P — P) is con- 
served. This property is nothing but the conservation of 
lepton number found in the standard case [H, |40|. In 
this situation, one expects to recover the standard dy- 
namics. In the following, we shall thus concentrate on 
the cases where g is not parallel to L, i.e., when flavor- 
changing currents are present. In the standard case, 
where B • D = ~Dz (for 9 = 0) is preserved, to get a 
diagnostic of the system it is enough to follow e.g. the 
evolution of Pz , which describes the polarization state of 
antineutrinos. In the present case, it is useful to mon- 
itor also Dz, where the differences with respect to the 
Standard Model case manifest more clearly. 



III. FLAVOR EVOLUTION 
A. Vacuum case 

To illustrate the effect of nonstandard neutrino- 
neutrino couplings on the flavor evolution of a neutrino 
gas we first consider a case in which the matter term 
A = 0. Since we have chosen 9 = 0, we assume that only 
51 ^ without loss of generality. We start our discussion 



of NSI effects referring to values of |g| < 10~^. In this 
regime, the terms Odgp) in Eq. ^ are sub-leading and 
it is clear that the new terms have the form g x Pp. 

The nonstandard couplings are expected to play a sig- 
nificant role when the transverse nonstandard term ^fJ-rQi 
in Eq. ^ dominates over the vacuum term loBz = —to. 
Thus, we expect NSI effects to be relevant for interac- 
tion strengths gi ^ 10~^. For smaller values, we expect 
to recover the standard picture. However, we note that 
also smaller values of nonstandard coupling can have an 
interesting consequence. If we assume vanishing vacuum 
mixing, we do not expect any flavor conversion in the 
standard case, since the polarization vectors are exactly 
aligned with B. However, in this situation, the nonstan- 
dard term transverse to B is enough to trigger an in- 
stability in flavor space, producing a small offset with 
respect to the B direction. This is enough to obtain fla- 
vor transformations in inverted hierarchy analogous to 
the standard case (39| . In particular, we initially observe 
the synchronized oscillations in which the polarization 
vectors essentially remain fixed to their initial value un- 
til the start of the bipolar oscillations, which lead to the 
inversion of the polarization vectors, conserving lepton 
flavor number at each step. In normal hierarchy the sys- 
tem is initially in its stable configuration and the non- 
standard effects have no impact. The behaviour of Pz is 
represented in the top left panel of Fig. [T] for the case of 
91 - 10-^. 

For 5i ^ 10 we observe that a transverse nonstan- 
dard component is enough to produce wild oscillations 
in the system. Essentially, the nonstandard interactions 
act as an external force on the standard system, violat- 
ing the symmetry that keeps Dz constant. Naively, it 
seems that the effect can be easily accounted for by go- 
ing into a frame co-rotating along g [iO,!!!]: we have just 
introduced a nonstandard matter potential profile. The 
angle between B — now rotating — and g is not necessarily 
small if g has a component perpendicular to B, even for 
small vacuum mixing angle. More important, however, 
is that in this case the initial condition of the neutrino 
state might be significantly offset compared with the g 
vector, i.e., P(0) and P(0) might have a significant com- 
ponent orthogonal to g, resulting in a large amplitude of 
oscillation. 

In presence of transverse nonstandard terms, the pos- 
itive electron lepton number we started with is not con- 
served, since 



bz ~ MrC [giDy - g2Dx\ 



(9) 



In Fig. [2] (top left panel) we show the behaviour of Dz 
for the case of gi — lO^'^. At small radii Dz oscillates 
with a full amplitude. However, reducing /i^ the ampli- 
tude of the oscillations diminishes: The decline of the 
upper envelope of Dz follows the decrease of gi^r com- 
pared to the vacuum oscillation frequency uj. Referring 
to the "pendulum" analogy [H, H^, the external force 
introduced by the NSI is initially strong enough to com- 
pletely dominate the dynamics. However, as the strength 
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FIG. 1: The single-angle evolution of Pz in normal (dashed curves) and inverted hierarchy (solid curves) for 9 — in vacuum 
(left panels) and in presence of matter (right panels). We consider values of nonstandard coupling gi equal to 10"^ (top panels), 
10~^ (central panels) and 1 (bottom panels). Note the different scale for r in the lower left panel. 



g^r of the external force decreases, it is finally so much 
weaker than the potential lu that it can no longer lead 
to significant oscillations. If the NSI term is sufficiently 
strong, the excess in electron neutrinos we started with 
will relax to complete flavor equilibration. 

In Fig. [3l we show the behaviour of at the end of 
the flavor conversions as a function of gi. We observe 
that nonstandard effects start to arise for gi > 10~*, in 
agreement with our prediction. For 10~^ ^ 5i ^ 10"^ 



we see that Dz saturates at zero. The corresponding be- 
haviour of Pz is shown in Fig. 1 for gi = 10~^ (central 
left panel) for both normal and inverted hierarchy. At 
the beginning, it oscillates with full amplitude which is 
then reduced following the decrease of fir- When the 
NSI effects have saturated and Dz has reached zero, the 
polarization vectors follow the standard "pendulum" dy- 
namics with null lepton number ^3§\. As /i^ — > 0, the 
kinetic energy of the flavor pendulum tends to zero and 
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FIG. 2: The evolution of Dz in normal (dashed curve) and inverted (solid curve) in vacuum (left panels) and in presence of 
matter (right panels). We use gi = 10^^ (upper panels) and g\ = 1 (bottom panels). Note the different scale for r in the lower 
left panel. 




log(g) 

FIG. 3: The value of at the end of flavor evolution as 
a function of g\ for inverted neutrino mass hierarchy in the 
vacuum case (solid curve) and in presence of matter (dashed 
curve). Similar to the normal mass hierarchy case (not 



in both normal and inverted hierarchy the polarization 
vectors tend to minimize the potential energy [39| 

£;. = |(|S|+B.S) , (10) 



where S = P + P. In normal neutrino mass hierarchy, 
where points in the negative z-direction, Pz rises 
again and returns to its initial value. In inverted hi- 
erarchy, where B points in the positive z-direction and 
the original equilibrium position was unstable, the ini- 
tial dynamics are enough to trigger significant flavor 
changes. The resulting behaviour is the total inversion 
of Pz- The main effects of a sub- leading nonstandard 
neutrino-neutrino interaction is thus to produce a flavor 
equilibration in the ensemble. 

We now consider large values of the g couplings. Even 
if from a theoretical perspective these are quite unnat- 
ural, the loose experimental bounds on these interac- 
tions allow us to play also with these wilder scenarios. 
Large couplings produce a modulation of the standard 
neutrino- neutrino term in Eq. ([6]). Note that in the ex- 
treme case of |g| = 2, the NSI term exactly compen- 
sates the standard interaction term. If it were not for the 
other nonstandard term on the right hand side of Eq. 
the system would become linear. This term produces a 
precession of the polarization vector P around the fixed 
g direction in flavor space, with a speed modulated by 
/irg-D. For generic large values of the NSI couphngs, the 
motion of the polarization vectors in flavor space is thus 
a combination of two precessions, one around D and the 
other around g, with comparable speed. Under this con- 
dition, the final value of at the end of the evolution is 
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of Pz and Dz as functions of r is analogous to the stan- 
dard case as shown in Figs. [T] and (upper right panels), 
respectively, for the case of gi — 10^^. If we allow for 
gi ~ 0{1), it can dominate over the A term close to the 
neutrinosphere. The behaviour of and Dz as functions 
of r for 51 = 1 is shown in Figs. [T] and [2] (lower right pan- 
els), respectively. We observe that, for gi = 1, the vectors 
initially undergo very rapid oscillations. However, when 
firdi < A, the standard matter potential dominates the 
dynamics and the oscillations only occur around the fla- 
vor direction, thus decreasing the flavor oscillation am- 
plitude. Thus, at the end of flavor evolution, the value 
of Dz saturates. 



FIG. 4: Radial profile for the matter potential A (solid curve) 
and for the neutrino- neutrino potential jj,r (dashed curve). 

quite unpredictable as shown in Fig. [3) It is determined 
by the detailed combination of the two different forces 
that drive the system: in particular, for gi — few 10~^, 
it is mostly the linear non-standard term that matters, 
while for larger couplings the quadratic corrections are 
important. The behavior of Pz and Dz as functions of r 
is shown in Figs. [1] and [5] (lower left panels), respectively, 
for the case oi gi = 1 in both normal and inverted hierar- 
chy. Initially, we observe wild oscillations in these vectors 
that decrease their frequency following the decline of fi^ ■ 
The evolution seems to maintain no track of the standard 
behaviour. The nice pendulum analogy, driving the po- 
larization vectors in different directions for the different 
neutrino mass hierarchies, remains. However, the fast 
oscillations of the final value of Dz as a function of gi 
inhibits the ability of making precise predictions. 



B. Effect of Matter Background 

In a realistic supernova environment, the propagation 
of a dense neutrino gas is also affected by the ordinary 
matter background. To characterize the matter potential 
A we refer to the time-dependent parametrization used in 
Ref. [H, [el] , inspired by supernova shock- wave simula- 
tions. In the following, we will focus on late times, fixing 
t = 7 s after the core bounce. In Fig. 21 we show the /i^ 
and A potentials for our input choices. 

For typical supernova conditions, the ordinary matter 
term AL dominates over sub-leading nonstandard inter- 
action terms. The effect of a dominant dense matter 
term is to project the EOMs along the weak-interaction 
direction [SSj. In the co-rotating frame around L, the 
fast-rotating g transverse component is enough to trig- 
ger the flavor conversions in inverted hierarchy in the case 
oi 9 — 0, but otherwise plays no crucial role. 

In particular, for our input choices, in Fig. [3] we ob- 
serve that the matter effect guarantees the conservation 
of Dz for gi ^ 0.3. For smaller values of gi, the behavior 



IV. CONCLUSIONS AND OUTLOOK 

In this work we have investigated the impact of non- 
standard neutrino-neutrino interactions on the flavor evo- 
lution of a dense neutrino gas. First, we have gen- 
eralized the equation of motion for a neutrino ensem- 
ble in presence of nonstandard neutrino-neutrino inter- 
action terms. Then, we have shown that if ordinary mat- 
ter effects are sub-dominant, flavor-changing neutrino- 
neutrino couplings lead to a dynamics very similar to 
the Standard Model analogue, with the additional fea- 
ture of a complete flavor equilibration of the neutrino 
ensemble, both in normal and in inverted hierarchy. In- 
cluding matter effects of the magnitude expected in a 
realistic supernova environment, sub- leading NSI correc- 
tions are suppressed by the dominant dense matter ef- 
fect and a picture similar to the "standard" one emerges. 
Small NSI thus have negligible effects, except for trigger- 
ing flavor conversions in inverted hierarchy even in the 
limiting case of 6* = 0. All in all, we conclude that the 
phenomenon of collective neutrino flavor evolution in su- 
pernovae is generally robust with respect to sub-leading 
neutrino-neutrino interactions of the four-fermion type. 
Potentially large effects can arise for values of the non- 
standard neutrino-neutrino couplings comparable to the 
ordinary ones, which are phenomenologically allowed but 
theoretically disfavored. 

Still, to make the problem tractable we made several 
approximations compared to a realistic physical system. 
Let us conclude this section by speculating on possible 
consequences of our findings in more realistic environ- 
ments. For example, in some situations the ordinary 
matter effect might be weak enough for sub-leading NSI 
effects to affect the dynamics. Perhaps the most notable 
example is in 0-Ne-Mg SN progenitors with rapidly de- 
creasing density profiles, where in particular the neutron- 
ization burst might be affected [H, [4^[65l|. Also, pecu- 
liar time-dependent signatures in the supernova neutrino 
burst could be produced by non-negligible NSI effects, 
since, at least at sufficiently late time after core bounce, 
the nonstandard neutrino-neutrino term could dominate 
over the dense matter term in the low-radii region. An- 
other physically interesting environment is provided by 
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accretion disks formed by coalescing neutron stars, where 
emitted neutrinos propagate in a region of relatively low 
matter density along the disk axis [s^l- Nonstandard 
interactions could play some role in this context, per- 
haps with an impact on the gamma-ray burst production 
by neutrino-neutrino annihilations in such engines. Also 
note that a more realistic treatment of the problem might 
require a detailed three-flavor study, which goes beyond 
our present purpose. 

Further directions of investigation may be related to 
the presence of nonstandard neutrino-charged fermion in- 
teractions. Their impact on matter transitions in super- 
novae has been explored in different works. In particular, 
in Ref. [13] it has been shown that these could induce new 
effects just above the neutrino sphere. Their impact on 
the collective neutrino flavor evolution still remains to be 
characterized. 

Finally, we would like to remark that this study as- 
sumes that the nonstandard effects do not affect the for- 
mation of the neutrino spectra and the dynamics in the 
supernova core. It has been argued in ^^^^ suffi- 

ciently large NSI of neutrinos with SN matter may lead 



to drastic changes in the core electron fraction and neu- 
trino spectra. Since the NSI considered here do not enter 
directly into the weak equilibrium, nor enhance neutrino- 
nucleus collisions, one would naively guess that at least 
for sufficiently small g their effects should be sub-leading. 
Of course, a definitive proof of this working hypothesis 
would require detailed supernova simulations, including 
nonstandard interaction effects in the dynamics of the 
core-collapse. 
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